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ABSTRACT 



We present deep composite luminosity functions in B, g,V , r, i and z for six 
clusters at 0.14 < z < 0.40 observed with the Hubble Space Telescope Advanced 
Camera for Surveys. The luminosity functions reach to absolute magnitude of 
~ — 14 + blogh mag. and are well fitted by a single Schechter function with 

MsgWiz = -19-8> -20.9 -21.9, -22.0, -21.7, -22.3 mag. and a 1.3 (in all 

bands). The observations suggest that the galaxy luminosity function is domi- 
nated by objects on the red sequence to at least 6 mags, below the L* point. 
Comparison with local data shows that the red sequence is well established at 
least at z ^ 0.3 down to ~ 1/600*'' of the luminosity of the Milky Way and that 
galaxies down to the regime of dwarf spheroidals have been completely assem- 
bled in clusters at this redshift. We do not detect a steepening of the luminosity 
function at M > —16 as is observed locally. If the faint end upturn is real, the 
steepening of the luminosity function must be due to a newly infalling population 
of faint dwarf galaxies. 

Subject headings: galaxies: luminosity function, mass function — galaxies: dwarfs — 
galaxies: formation and evolution — galaxies: clusters 
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Introduction 



The galaxy luminosity function (hereafter LF) provides a 0*'^ order description of 
the gross properties of galaxy populations. Although it may be regarded as somewhat 
simplistic, in reduci n g gala xy properties to the two numbers (M* and a) that describe the 



LF in the 



Schechten (119761 ) form, reproducing the observed LF is a fundamental test for any 



viable theory of galaxy forma tion, and one that has proven surprisingly complex to achieve 



until relatively recently (e.g., 



Croton et al. 



20061). 



Clusters of galaxies are, in several ways, ideal laboratories to study the LF and its 
evolution. Cluster galaxies may be considered as a volume limited sample of objects 
observed at the same cosmic epoch and occupy a relatively 'constant' environment 
corresponding to the highest density peaks at the redshift of observation. Operationally, 
cluster members have a much higher surface density on the sky than the surrounding field 
and often exhibit distinctive morphologies and colors (see Fig. 1 below). It is therefore 
possible to establish membership in a cluster (at least in a statistical sense) without 
resorting to observationally expensive redshift surveys. 

On the other hand, clusters are obviously special places, containing only a small 
percentage (5-10%) of all galaxies in the Universe. Their galaxy populations clearly 
implicate the cluster envir onment in a va riety of processes that affect the morphological 



evolutio n of galaxies (e.g. 



history (ILewis et al 



2002 



Dressier 



Gomez et al 



19801) an d suppress or modify their star formation 



20031 ). In order to use clusters as probes of galaxy 



evolution it is therefore important to understand the effects of environment and carry out 



careful comparisons with field sam ples. For loca . 



large redshift surveys such as 2dF (IColless et al. 



sam ples, this is now possib. 



200ll) and SDSS flYork et al. 



e thanks to 



2000). 



Although bright galaxies appear to have forme d the majority o: 
and assembled their mass at least by 2; = 1.5 (e.g. 



Andreon 



2006a 



^heir stellar popu 



De Propris et al. 



ations 



2007 : 
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Muzzin et al.ll2008l and references therein), there is considerable evidence that fainter 
galaxies undergo significant evolution since z = 1. The red s equence in clusters in the 



EDisCS sample ap pears to be truncated at faint magnitudes (jPe Lucia et al. 



Stott et al. 



20071 ). while 



(120071 ) sho w that the 



higher redshifts, as do 



Krick et al 



uminosity function of red cluster galaxies flattens at 



(120081) for clusters in a deep Spitzer field. On the other 



hand, there are c 



(lAndreon 



2006b 



usters with we 



Crawford et al 



formed red sequences at redshifts approaching 1 as well 



20081 ). It may not be surprising that the faint end of the 
LF shows considerable variations from cluster to cluster, as one expects that dwarfs are 
more sensitive to environmental effects. Although the bright end of the LF (M < —16) 
is broadly univ ersal, the faint end yaries between clusters and may even change radially 



within clusters (jPopesso et al. 



2006 



Barkhouse et al. 



2003). 



The evolution of the faint end slope of the total LF provides an interesting test of 
galaxy formation models. Even if the red sequence is truncated at faint luminosities 



at earlier epochs, one expects a steepening fai nt end as estab 



spectrum of fluctuations, which is very steep (IKhochfar et al. 



ishec 



by the initial power 



20071 ). In higher redshift 



clusters, there should be a steeper LF and a larger fraction of bluer dwarfs. 

The purpose of this paper is to explore the evolution of the LF to very faint magnitudes 
in a sample of z ~ 0.3 clusters observed with the Hubble Space Telescope (HST), reaching 
well below the L* point and into the domain of dwarf spheroidal galaxies, using a multicolor 
imaging sample. This will allow to test scenarios for the evolution of dwarf galaxy 
populations in ACDM hierarchical models. The structure of this paper is as follows: we 
describe the data and photometry in the next section. We derive the composite LFs 
and then discuss the results in the context of galaxy formation models. We adopt the 
'consensus' cosmology with Qm = 0.27, Qa = 0.73 and use Ho=100 km s^^ Mpc^^. All data 
are photometrically calibrated to the AB system, using the latest ACS zeropoints published 
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on the HST Instrument Web page. 



2. Observations and Data Analysis 



Observations for this paper consist of deep HST imaging of six clusters at 
0.14 < z < 0.40 taken with the Advanced Camera for Surveys (ACS) in the B (F435W), g 
(F475W), V (F555W), r (F625W), i (F775W) and z (F850LP) bands with exposure times 
of 5,000-10,000 seconds in each band. Each cluster was observed in a single ACS exposure, 
covering about 200" on the sky. The data were originally taken for studies of gravitational 
lenses in these clusters. Table 1 shows the clusters observed, their redshifts, exposure 
times in each band, and HST proposal IDs. Table 2 summarizes some essential physical 



charac teristics for the observed clusters. Most data are fron a the compilatio r i of 



( 119991 ). from which r2oo can be estin aated via the formula of 



for A1413 where r2oo was derived by 



Pointecouteau et al. 



Carlberg et al. 



Wu et al. 



(119971 ). except 



(120051 ) from the X-ray profile. We 



were unable to find data in the literature for A1703. 



The data were retrieved from the HST archive as fully red uced and fiatfielded files 



2OO2I) to 



(*.fit) and then processed through the Multidrizzle algorithm (iKoekemoer et al.l 
produce fully registered and co-added images for each band. Fig. 1 shows color composites 
for our data. Full resolution JPEGs will be made available on the AJ web site. 



Detection of objects and photometry were carried out using Sextractor (IBertin fc Arnouts 



19961 ). We experimented with Sextractor parameters to maximize our detections and 
minimize the number of spurious objects. All detections were visually examined to eliminate 
noise spikes, bleed trails from bright stars and other contaminants, especially the numerous 
arclets present in the images. The Sextractor parameters eventually employed are shown in 
Table 3. We used the same parameters for all images. 
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For each object we compute a total (Kron-like) magnitude and an aperture magnitude 
of area equivalent to the minimum number of connected pixels needed for detection. This 
provides an estimate of the central surface brightness for each object. The motivation 
behind this procedure is as follows. Detections of objects in deep images depends on their 
total magnitude but also on their central surface brightness. An object could be brighter 
than the magnitu de threshold but be lo st in the night sky because of low surface brightness 



fsee discussion in 



Cross fc Driver 



2OO2I ) . For this reason we need to determine both a 



limiting total magnitude for completeness and a surface brightness threshold for detection. 

We do this by plotting central surface brightness vs. total magnitude for all bands 
in Fig. 2. To save space we only show the data for Abell 1703, but equivalent figures for 
all other clusters are available on the AJ web site. In these figures the sequence at high 
surface brightness that provides an upper limit to the scatter plot is ca used by stars. We 



Garilli et al. 



I999I). 



can therefore discriminate against stellar contamination in this way (cf. 
The surface brightness threshold is chosen empirically. The limit is selected to include as 
many objects as possible, but avoiding regions of the parameter space where the sample is 
obviously very incomplete. Selection lines in surface brightness, apparent magnitude and 
the star- galaxy separation line are shown in Fig. 2. 

In order to determine a limiting magnitude, we plot the raw number counts for 
galaxies in the Abell 1703 field in Fig. 3. As for Fig. 2, figures for all the other clusters 
are made available electronically. The completeness magnitude is chosen to lie about 0.5 
mag. brighter than the luminosity at which counts start to decrease, in order to select a 
highly complete sample of objects. 

The only way to establish cluster membership for faint galaxies is by statistical 
background subtraction. In order to do so, we need to analyze 'blank' (cluster-less) 
comparison fields of similar or greater photometric depth. We choose to use the two 
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GOODS fields (iGiavalisco et al.ll2004l ) as these are the deepest and widest fields where HST 
multicolor photometry is available. 

We used the same Sextractor parameters as we used for the cluster fields and impose 
the same selection limits on the GOODS fields as we did for the target fields. Fig. 4 shows 
the detection and selection plots for the GOODS fields (similar to Fig. 2 above). The fields 
are much deeper than the cluster fields we use. It must be noted that we plot only 10% of 
the detections in the GOODS fields in Fig. 4 to avoid saturating the figure. 

Fig. 5 plots galaxy number counts for the GOODS fields and comparison galaxy number 



counts from the literature. The HST data were corrected to the SDSS and 



systems using the transformations tabulated in IHolberg fc Bergeron! (120061 ). GOODS 



ohnson- Cousins 



number counts were fitted with a quadratic of the form oq + aix + a2X^. The fit was carried 
out over galaxies brighter than the completeness limit shown in Fig. 4 and fainter than 
m~ 21 mag. as GOODS counts for brighter galaxies are affected by small scale clustering. 
Table 4 shows the values of the coefficients for the quadratic fits. 



3. Luminosity Functions 

We now subtract the scaled contribution from fore/background counts from galaxy 
counts. Errors in galaxy counts for clusters and the GOODS fields are assumed to be 
Poissonian. We used our smoothed quadratic for the background counts, extrapolating for 
galaxies brighter than m ~ 21 mag. Field counts at these bright limits and over the small 
(~ 9 arcmin^) covered by ACS are very small in any case. Clustering errors over the cluster 
and GOODS fields a re estimated f ollowi ng the counts-in-cells approach and the fitting 



formula described by iHuang et al 



fll997h . 



It may be possible to further refine the selection of cluster members using photometric 
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redshifts or color cuts. Contamination of the LF by backgr ound galaxies has be en suggested 



as a possible cause of the faint end upturn (see below) 
this would lead to loss of gen erality. The argumen t by 



Valotto et al 



Valotto et al. 



fIgOOlh . However, 



for nearby clusters (although 



Andreon et al. 



2005 



(I2OOII ) may be correct 



Valotto et al. 



2001 



point out that the simulations used by 



are unrealistic and that the background subtraction approach used here 
works as long as proper care is taken of statistical uncertainties, see their §3.3.6), but it is 
not appropriate for these distant objects where the redshift is much larger than the typical 
scale of structure in the Universe (so that the counts should be smooth). 

There are some mismatches between the GOODS data (used for background 
subtraction) and the cluster data. For the g (F435W) and r (F625W) bands, we can use the 
GOODS F475W and F606W counts (respectively), which are close matches to the bands 
used, however, for F555W {V) we used the average of the F475W and F606W counts. The 
LF in this bands is more uncertain, because we do not have an appropriate background 
field. 

We k and e c orrect the LFs to z = fo r ease of comparison with local data. In order 
to do so we use a 



Bruzual &: Chariot 



(120031 ) model for a solar metallicity elliptical formed 
at 2; = 3 with an e-folding time of 1 Gyr. This is appropriate to the red ellipticals that 
dominate the bright end of cluster LFs. In effect this approach allows us to consider the 
differential evolution, if any, between the dwarfs and the giants, who appear to consist 
mainly of old galaxies with little or no recent star formation. 

The LF of each individual cluster is rel atively u ncert ain. We therefore build composite 



LFs in each band following the method by 



CollessI (119891 ) in order to average out errors. 



The method assumes that the LF is broadly universal, but this appe ars to be borne out 



by observations at low redshift, at least for relatively bright galaxies (jPopesso et al. 



2006 



Barkhouse et al. 



20071 ). 



-9- 



Composite LFs are built as follows: the number of cluster galaxies in the j*'* magnitude 
bin of the composite LF is given by: 



where Nij is the (background corrected) number of galaxies in the j*^ magnitude bin of 
the LF of the i^^ cluster LF, N^q is the normalization of the cluster LF (to take care of the 
different richnesses of the clusters) and is taken to be the background corrected number of 
cluster galaxies brighter than Mb = — 19 (or its equivalent in the other bands, assuming the 
color of an old elliptical), ruj is the number of i clusters contributing to the magnitude 
bin of the composite LF and Ncq is the sum of all the i normalizations. 



Essentially, this carried out a weighted average of the individual cluster LFs, scaled 
by the number of galaxies (in each cluster) brighter than approximately the L* point. The 
formal errors on N^j are computed as: 



where 5Nij are the errors on the number counts (including Poissonian and clustering 
contributions) for the magnitude bin of the i*^ cluster LF. 

For each cluster we build the composite LF by summing up the individual LFs in 
absolute magnitude bins, after correcting for the distance modulus, extinction and k -\- e 
correction to the common redshift of 2; = 0. 

The resulting LFs are presented in Fig. 6, together with the 1, 2 and 3(7 error countours. 




(1) 




(2) 




(3) 
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We fit the data to a single Schechter function (see discussion below, regarding the existence 
of a faint end upturn). In some cases we exclude the fainter and/or brighter magnitude bin 
from the fit, because these points are dominated by either the poorly sampled very massive 
galaxies or the fainter objects that suffer from complex incompleteness due to magnitude 
and surface brightness selection. The fit parameters are tabulated in Table 5, together with 
the la errors for each parameter, marginalized over the remaining parameters. 



4. Discussion 



There are only a few deep cluster LFs having comparable depth to our own. In general. 



(Driver et al. 


1994 


r 


Crent 


lam 


1998a: 


. — ^ — „^ 
Bovce et al. 


2001; 


2004; 


Andreon et al. 


2005 




Harsono & De Propris 


2007) 



Mer curio et al 



lallower samples 



2003 



Pracy et al. 



20071 ). The M* point we measure gets 



steadily brighter with increasing wavelength, tracking the color of old, passiv ely evolving 



stella r populations that are the main component of bright elliptical galax ies (ICool et al. 



20061) . The value we find is also consistent with the local SDSS value by 



Po pesso et al. 



( 120061 ): this is not surprising, as it is well known that massive galaxies evolve passively since 



high redshift. 



The main interest of this investigation is in the evolution of the LF slope. The LF is 
well fitted by a single Schechter function, with a slope a ~ —1.3 which is essentially the 
same in all bands. This suggests that the cluster LF is dominated by galaxies on the red 
sequence to a luminosity of M ~ —14 in all bands, and that therefore the red sequence 
continues at least to 6 magnitudes below the L* point even at z ~ 0.3. This can be clearly 
seen in Fig. 7 where we show the g — r vs. z color-magnitude relation for Abell 1703. The 
relation can be clearly followed at least to 2; = 25 with little scatter. A more detailed 
analysis of the color-magnitude relation is deferred to a future paper (Harsono & De Propris 
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2008, in preparation). The red sequence appears to be well established an d contain the 



majority of cluster populatio ns even at a lookback time of ~ 3.5 Gyrs (cf. 



Andreon et al. 



2006 



Eisenhardt et al. 



20071 for a similarly deep relation in local clusters). 



he sl ope is very simila r to th e local value for M < — 16 derived by 



(120061) and 



Barkhouse et al. 



Popesso et al. 



(120071 ). This argues that even galaxies to M ~ —14 have 



formed their stellar populations and assembled the majority of their mass at least at 
z ~ 0.3. We can therefore place a significant lower limit to the assembly epoch of galaxies 
down to 1/600*^ of the mass of the Milky Way. Given the evidence for truncatio n of the 



red sequence at z ~ 0.8 (jPe Lucia et al 



2007 



Stott et al 



2007 



Krick et al. 



20081 ) ■ our data 



point to the 0.3 < z < 0.8 interval as a crucial epoch to investigate the formation of the 
fainter galaxy populations. 

The other issue we wish to address is t he faint end upturn . Thi s has been a contr o versia l 



subject, ever since its original discovery by 



Driver et al 



fll994h and 



De Propris et al 



(119951). 



There have been a series of claims and counterclaims regarding the faint end slope and 
the existence of an upw a rd inflection, some ti mes even in the same cluster. For instance. 



De Propris et al. 



and 



Milne et a. 



Bernstein et al 



19981); 



(119951); 



TrenthamI fll998bh 



Adami et al. 



(12007h : 



Jenkins et al. 



(120071) 



(l2007h observe a ste e p LF at t he faint end for the Coma cluster, while 



Adami et al 



flatter slopes. I n Vir go 



Rines fc Gellei] (120081 ) derive a = -1.3 



Sabatini et al. 



(Il998h a nd 



Beijersbergen et al. 



(I2OO2I ) quote much 
(120051) claim a slope as steep as ~ —1.6, while 



Baldry et al. 



(I2OO8I ) review the existence of the 



upturn and the shape of the faint end of the LF and conclude that there is a steep mass 



function. The observation o: 



groups (jPopesso et al. 



2006 



the upturn in two co mposite LFs, derived by two different 



Barkhouse et al 



the LF is real (however, cf., iHilker et al 



2003 



20071). suggests that the u pward inflection of 



Penny fc Conselice 



20081 ). 



Unlike 



Popesso et al. 



(120061 ) and 



Barkhouse et al. 



(I2OO7I ) we do not find a steepening 
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of the LF at M > —16. However, we must consider that the LF upturn is more pronounced 
at large clustercentric radii, while our fields generally cover only the central region of each 
cluster. The size of the ACS fields used in this st udy covers between 350 and 750 kpc 



on the side. In terms of -R200 (the radius at which 



Popesso et al. 



2006 



and 



Barkhouse et al. 



(I2OO7I ) normalize their LFs), the areas imaged by ACS cover between 20% and 40% of the 
area of the cluster out to r2oo- We therefore only derive LFs for the cluster cores. This 
should not affect our comparisons with the bright end of the LF (M < —16), a s this does not 



seem to va ry significantly wi t 



(120061) and 



Barkhouse et al. 



1 rad ius. However, the steep upturn claimed by 



Popesso et al. 



(120071 ) is much more pronounced at large clustercentric radii. 



i.e., in the cluster outskirts. On the other han d it is possible to see the upturn even in the 



more cent ral regions of t 



Fig. 10 of 



Popesso et al. 



le sa mple studied by 



Popesso et al. 



(l2006l ). inside ~ 0.3r2oo (see 



20061 ). so our data should exhibit an upturn at faint magnitudes. 



although of course more clusters would be welcome to bolster the argument. 

As we do not detect a steepening in the LF, we suggest that, if the steepening is 
real, the faint galaxie s contributing to the upturn consist of a population of recently 



infalling objects (e.g.. 



Smith et al. 



20081 ) whose star formation is curtailed by the cluster 
environment. This would be consistent with the observation that about 1/2 of the fainter 
cluster dwarfs in Virgo and elsewhere have been forming stars until relatively recently 



(IConselice et al. 



2001 



2OO3I ). until their star formation stopped and their colors moved 



to the red sequence. The newly infalling population may be the one now contributing to 
the steep upturn. In a 'downsizing' model, it may be expected that objects undergoing 
star formation suppression and cluster infall at the present epoch would indeed tend to 
be among the fainter and less massive galaxies and to have the steep LF characteristic of 
'pristine' CDM power spectra. 



Our observations imply that the passive evolution of bright galaxies can be extended 
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to faint dwarfs, at least to 2; ~ 0.3 and suggest that the majority of galaxy evolution may 
have taken place at surprisingly early epochs even for the least massive objects. 

We would like to thank the referee for a very helpful report, that has considerably 
improved the paper. We would also like to thank Stefano Andreon for having commented 
on a version of this paper. 
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at the Space Telescope Science Institute (MAST). STScI is operated by the Association of 
Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. Support 
for MAST for non-HST data is provided by the NASA Office of Space Science via grant 
NAG5-7584 and by other grants and contracts. 



We acknowledge use of the Cosmology Calculator (IWrightl 120061 ). 
Facilities: HST (ACS). 



- 14- 
REFERENCES 

Adami, C, Nichol, R. C, Mazure, A., Durret, F., Holden, B. & Lobo, C. 1998, A&A, 334, 
765 

Adami, C, Picat, J. P., Durret, F., Mazure, A., Pello, R. & West, M. 2007, A&A, 472, 749 
Andreon, S., Punzi, G., Grado, A. 2005, MNRAS, 360, 727 
Andreon, S. 2006a, A&A, 448, 447 
Andreon, S. 2006b, MNRAS, 369, 969 

Andreon, S., Cuillandre, J.-C, Puddu, E. & Mellier, Y. 2006, MNRAS, 372, 60 

Baldry, I. K., Glazebrook, K. & Driver, S. P. 2008, MNRAS, 388, 945 

Barkhouse, W. A., Yee, H. K. C. & Lopez-Cruz O. 2007, ApJ, 671, 1471 

Beijersbergen, M., Hoekstra, H., van Dokkum, P. G.& van der Hulst, T. 2002, MNRAS, 
329, 385 

Bernstein, G. M., Nichol, R. C., Tyson, J. A., Ulmer, M. P. & Wittman, D. 1995, AJ, 110, 
1507 

Bertin E., Arnouts S. 1996, A&AS, 117, 393 

Boyce, P. J., Phillipps, S., Jones, J. B., Driver, S. P., Smith, R. M., & Couch, W. J. 2001, 
MNRAS, 328, 277 

Bruzual, G. & Chariot, S. 2003, MNRAS, 344, 1000 

Capak, P. et al. 2007, ApJS, 172, 99 

Carlberg, R. G., Yee, H. K. C. & EUingson, E. 1997, ApJ, 478, 462 



- 15 - 

Crawford, S. M., Bershady, J. A. & Hoessel, J. G. 2008, astro.ph 0809.1661 
CoUess, M. M. 1989, MNRAS, 237, 799 
CoUess, M. M. et al. 2001, MNRAS, 328, 1039 

Conselice, C. J., Gallagher, J. S. & Wyse, R. F. G. 2001, ApJ, 559. 791 
Conselice, C. J., Gallagher, J. S. & Wyse, R. F. G. 2003, AJ, 125, 66 

Cool, R. J., Eisenstein, D. J., Scranton, R., Brinkmann, J., Schneider, D. R & Zehavi, I. 

2006, AJ, 131, 736 

Cross, N. J. G. & Driver, S. P. 2002, MNRAS, 329, 579 
Croton, D. et al. 2006, MNRAS, 365, 11 
De Lucia, G. et al. 2007, MNRAS, 374, 809 

De Propris, R., Pritchet, C. J., Hartwick, F. D. A. & McClure, R. D. 1995, ApJ, 450, 534 

De Propris, R., Eisenhardt, P. R., Stanford, S. A. & Dickinson, M. 1998, ApJ, 503, L45 

De Propris, R., Stanford, S. A., Eisenhardt, P. R., Holden, B. P. & Rosati, P. 2007, AJ, 
133, 2207 

Dressier, A. 1980, ApJ, 236, 351 

Driver, S. P., Phillipps, S., Davies, J. I., Morgan, I., Disney, M. J. 1994, MNRAS, 268, 393 

Eisenhardt, P. R., De Propris, R., Gonzalez, A. H., Stanford, S. A., Wang, M. & Dickinson, 
M. 2007, ApJS, 169, 225 

GariUi, B., Maccagni, D. & Andreon, S. 1999, A&A, 353, 479 

Giavalisco, M. et al. 2004, ApJ, 600, L93 



-16- 



Gomez, P. L. et al. 2003, ApJ, 584, 210 
Harsono D. & De Propris R. 2007, MNRAS, 380, 1036 
Hilker, M., Mieske, S. & Infante, L. 2003, A&A, 397, L9 
Holberg, J. B. & Bergeron P. 2006, AJ, 132, 1221 

Huang J.-S., Cowie L. L., Gardner J. P., Hu E. M., Songaila A., Wainscoat R. J. 1997, ApJ, 
476, 12 

Jenkins, L. P., Hornschemeier, A. E., Mobasher, B., Alexander, D. M. & Bauer, F. E. 2007, 
ApJ, 666, 846 

Koekemoer, A., Fruchter, A. S., Hook, R. N. & Hack, W. 2002, in The 2002 HST Calibration 
Workshop: Hubble after the installation of ACS and the NICMOS cooling system 
ed. by S. Arribas, A. Koekemoer, B. Whitmore (Baltimore: Space Telescope Science 
Institute), p. 337 

Khochfar, S., Silk, J., Windhorst, R. A. & Ryan, R. E. 2007, ApJ, 668, L115 

Krick, J. E., Surace, J. A., Thompson, D., Ashby, M. L. N., Hora, J. L., Gorjian, V. & Yan, 
L. 2008, astro-ph 0807.1565 

Lewis, I. J. et al. 2002, MNRAS, 334, 673 

MacDonald, E. C. et al. 2004, MNRAS, 352, 1255 

Mercurio, A., Massarotti, M., Merluzzi, P., Girardi, M., La Barbera, F. & Busarello, G. 
2003, MNRAS, 408, 57 

Milne, M. L., Pritchet, C. J., Poole, G. B., Gwyn, S. D. J., Kavelaars, J. J., Harris, W. E. 
& Hanes, D. A. 2007, AJ, 133, 177 



- 17- 

Muzzin, A., Wilson, G., Lacy, M., Yee, H. K. C. & Stanford, S. A. 2008, astro-ph 0807.0227 

Penny, S. J. & Conselice, C. J. 2008, MNRAS, 383, 247 

Pointecouteau, E., Arnaud, M. & Pratt, G. W. 2005, A&A, 435, 1 

Popesso, P., Biviano, A., Bohringer, H. & Romaniello, M. 2006, A&A, 445, 29 

Pracy, M., De Propris, R., Driver, S. P., Couch, W. J. & Nulsen, P. E. J. 2004, MNRAS, 
352, 1135 

Rincs, K. & Geller, M. J. 2008, AJ, 135, 1837 

Sabatini, S., Davies, J. I., Scaramella, R., Smith, R., Baes, M., Linder, S. M., Roberts, S. & 
Testa, M. 2005, MNRAS, 357, 819 

Schechter, P. 1976, ApJ, 203, 297 

Smith, R. J. et al. 2008, MNRAS, 386, L96 

Stott, J. P., Smail, I., Edge, A. C., Ebehng, H., Smith, G. P., Kneib, J.-P. & Pimbblet, K. 
A. 2007, ApJ, 661, 95 

Trentham, N. 1998a, MNRAS, 295, 360 

Trentham, N. 1998b, MNRAS, 293, 71 

Valotto, C., Moore, B. & Lambas, D. G. 2001, ApJ, 546, 157 

Wright, E. L. 2006, PASP, 118, 1711 

Wu, X.-R, Xue, Y.-J. & Fang, L.-Z. 1999, ApJ, 524, 22 

Yasuda, N. et aL 2001, AJ, 122, 1104 

York, D. G. et al. 2000, AJ, 120, 1579 



-18- 



This manuscript was prepared with the AAS macros v5.2. 



-19- 



Fig. 1. — Colour images of target clusters. From left to right: A1413, A1689, A2218, 
A1703, MS1358 and C10024. These figures exceed arXiv size limits and can- 
not BE COMPRESSED FURTHER. THEY WILL APPEAR AS JPEG IN AJ. PLEASE EMAIL 
rdepropris@ctio.noao.edu FOR THE FIGURES TO BE SENT TO YOU VIA FTP (PROVIDE 
site) or email 
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Fig. 2. — Mean central surface brightness vs. total magnitudes for objects in the Abell 1703 
field. Selection lines show the star-galaxy separation, the surface brightness threshold and 
the magnitude completeness limit. Similar figures for all the other objects in our sample can 
be found in the electronic version of this paper. 
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Fig. 3. — Raw galaxy number counts in the Abell 1703 field. See the electronic version 
this paper for the relative figures for all other targets. 
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Fig. 4. — Same as Fig. 2, but for the combined GOODS fields. We only plot a random 
selection of 10% of the objects in the GOODS fields, for clarity. 
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Fig. 5. — Number counts for the GOODS fields and quadratic fits. Literature counts (see 
text) are also shown. 
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Fig. 6. — Composite luminosity functions, best fits and error ellipses for the B and g bands 
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, — (Continued): for the V and r bands 
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Fig. 6. — (Continued): for the i and z bands 
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Fig. 7. — The color magnitude diagram for galaxies in A1703 



Table 1. Archival HST cluster observations 



Cluster 


redshift 


B Exposiire 


9 


V 


r 


i 


z 


Proposal ID 


Abell 1413 


0.143 










2440 


2560 


9292 


AbcU 2218 


0.171 


7048 


5640 


7048 


8386 


10732 


5640 


9452, 9717, 10325 


Abell 1689 


0.183 




9500 




9500 


11800 


16600 


9289 


Abell 1703 


0.258 


7050 


5564 


5564 


9834 


11128 


17800 


10325 


MS 1358.4+6245 


0.328 


7928 


5470 


5482 


9196 


10952 


15017 


9292, 9717, 10325 


CI 0024.0+1652 


0.395 


6435 


5072 


5072 


8971 


10144 


16328 


10325 



Table 2. Cluster properties 



Cluster 


redshift 


(T [km 


s-i] 


T (KcV) 


Lx (10*4 [ergs s"!] 


r200 (Mpc /i-i) 


Coverage (% r2oo) 


Abell 1413 


0.143 






8.85 ± 0.50 


36.01 + 4.54 


1.71 


21 


Abell 2218 


0.171 


1800 ± 207 


9 02+" 


55.73 + 8.92 


2.94 


14 


Abell 1689 


0.183 


1370^ 


-180 
-220 


7.10 + 0.20 


21.96 + 2.01 


2.24 


19 


Abell 1703 


0.258 














MS 1358.4+6245 


0.328 


937 ± 54 


7.50 + 4.30 


21.81 + 3.81 


2.40 


28 


CI 0024.0+1652 


0.395 


1339 ± 233 




4.26 + 0.01 


1.95 


39 
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Table 3. Sextr actor parameters 



Parameter 


Value 


DETECT MINAREA 


7.0 pixels 


DETECT THRESH 


2.0 


ANALYSIS THRESH 


2.0 


DEBLEND NTHRESH 


32 


PHOT APERTURES 


1.5 


PHOT AUTOPARAMS 


2.5, 3.5 
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Table 4. Quadratic fits to GOODS counts 



Band 


aO 


al 


a2 


g (F475W) 


-44.46 


3.544 


-0.06346 


r (F606W) 


-23.58 


1.928 


-0.03196 


i (F775W) 


-18.51 


1.590 


-0.02643 


z (F850LP) 


-10.50 


0.9377 


-0.01417 
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Table 5. Best fit parameters for composite LFs 



Band 


M* 


a 




B (F435W) 


-19.80 ±0.27 


-1.29 ±0.06 


0.59 


g (F475W) 


-20.94 ±0.17 


-1.31 ±0.04 


0.95 


V (F555W) 


-21.86 ±0.27 


-1.27 ±0.05 


1.08 


r (F625W) 


-21.95 ±0.29 


-1.33 ±0.03 


0.43 


i (F775W) 


-21.66 ±0.27 


-1.27 ±0.04 


0.37 


z (F850LP) 


-22.26 ±0.30 


-1.45 ±0.02 


0.94 



